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ANALYTICAL COMPARISON OF EFFECTS OF SOLID -FRICTION 
AND VISCOUS STRUCTURAL DAMPING ON PANEL FLUTTER 

Herbert J. Cunningham 
Langley Research Center 

SUMMARY 

A Galerkin modal analysis is presented that accounts for the effects of both solid- 
friction and viscous structural damping on panel flutter, based on unsteady aerodynamic 
forces from supersonic potential flow. The eigensolutions are made by complex eigen- 
value computer routines. 

Markedly different effects on the flutter boundary of the two types of structural 
damping are obtained. This result establishes that there is not, in general, an "equivalent 
viscous" damping for solid-friction damping. 

For the limiting case of the static-aerodynamic approximation, a substantially dif- 
ferent flutter dynamic pressure is obtained for solid friction identically zero compared 
with solid friction approaching zero as a limit. Use of the quasi-static aerodynamic 
approximation eliminates that difference. 

INTRODUCTION 

Two representations of structural damping have been used in published panel-flutter 
analyses. The two types or representations are: (1) solid friction that is termed 

"hysteretic" in some of the literature and is accounted for within a complex stiffness; 
and (2) the viscous type, that for harmonic motion is also characterized by a hysteresis 
loop in stress-strain and load -deflection plots. 

Neumark in reference 1 has provided a thorough and careful analysis of the basic 
characteristics of viscous and solid-friction (hysteretic) damping for single -degree -of - 
freedom systems in free and forced oscillation; his findings apply as well to multiple- 
degree-of -freedom systems oscillating in a single simple harmonic motion. Perhaps the 
earliest application of viscous-type structural damping to panel-flutter analysis appears 
in reference 2. The model analyzed was an infinite streamwise array of panels with a 
continuous spectrum of travelii^-wave modes, and some ranges of parameters were 
found for which added dampii^ was destabilizing. 

Reference 3 included an early application to panel -flutter analysis of the solid- 
friction type of structural damping used in wing -flutter analysis in reference 4. Added 


damping was found to be stabilizing in some parameter areas but destabilizing in others. 
An analysis is given in reference 5 of plane panels subjected to in-plane loading, linear 
piston-theory aerodynamics, and with solid-friction (hysteresis) damping both in bending 
and in longitudinal compression. The destabilizing effect of structural damping for some 
ranges of parameters was found also for this case. Reference 5 is notable for its com- 
prehensive list of references (69). 

The extensive investigation of reference 6 included some consideration of viscous- 
type structural damping and found that it had a moderate destabilizing effect for some 
parameters, especially in the presence of low aerodynamic damping that accompanies low 
air density. A careful study is given in reference 7 of the effects of both viscous and 
hysteretic (solid -friction) damping on the flutter of a rectangular membrane. It was 
found that an increase of either or both types of damping is destabilizing, and that, for the 
same value of structural damping coefficient, viscous damping is more destabilizing than 
solid-friction (hysteretic) damping. 

The purpose of the present report is to present two methods of solution for routine 
calculation of panel-flutter characteristics as affected by the two types of damping in the 
presence of the complete unsteady aerodynamic forces and to show the substantially con- 
trasting effects of the two dampings. The present work is an outgrowth from that reported 
in reference 8 for unyawed panels. A yawed panel is analyzed herein to show the effects 
from the two dampings. In general, the complete unsteady aerodynamic forces from 
potential -flow theory are employed. The static and quasi-static aerodynamic approxima- 
tions are also used for comparison purposes. The equations for the generalized aero- 
dynamic forces are given in an appendix. 

SYMBOLS 

Et^ 

D flexural stiffness of panel, — ■, where v is Poisson’s ratio and 

12(1 - v^) 

E is Young's modulus 

ggj modal -independent coefficient of solid -friction (sf) structural damping 

(eq. (12c)) 

gi g^,gj coefficients of solid-friction (sf) structural damping for modes i and j, 
respectively 

^j,sf ~^j,sf'^^sf 
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before being specialized to. mode j 


r 

®sf ^j,sf 

Si v’Si „}gr> ,, coefficients of viscous (v) structural damping for modes i, j, and B, 

respectively, where B denotes any base or reference mode 

Gg viscous structural damping constant (eq. (1)) 

hj = hj(x,y) mode-shape deflection for mode j 

Hj(x,y, t) mode -j component in deflection z (eq. (5)) 

i unit of imaginaries, fT 

k reduced frequency based on reference length £, cof/V 

Kj J nondimensional eigenvalue quantity for fundamental panel mode (see eq. (19)) 

£ panel length at zero yaw 

m^ mass of panel per unit surface area 

M Mach number of undisturbed stream 

Mj generalized mass for mode i (eq. (12)) 

Mq value of for h^ = 1.0; for \miform panels, Mq = m^£w 

p,q number of half-waves lengthwise and widthwise, respectively, in panel 

modes hj 

Ap perturbation pressure difference on panel, positive with z 

Apj complex amplitude of Ap for mode j (see eq. (8)) 

qj(T),qj time-varying generalized coordinate of motion for mode j and its complex 

amplitude, respectively, q.j(T) = e^^’’q. 

J J 

Qij,Qij generalized aerodynamic -force elements, Qy = Qjj^pV^fw^ (eq. (11)) 
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S surface area of panel, region of integration 

t panel thickness 

V speed of undisturbed airstream 

w panel width at zero yaw 

x,y lengthwise and widthwise coordinates, based on £ and w, respectively, 

fixed to the panel as it yaws 

z z-coordinates of deflected panel svirface 

^ - 1 


X 

A 

M 

P 

T 

0) 


CL)..Ct)wCO-r^ 

1’ B 


dynamic-pressure parameter, 
angle of yaw of panel, deg 


ratio of panel mass to air mass, 


mA^w _ Mq 
p£^w p£^W 



density of undisturbed stream 
time 

circular frequency of motion 

natural frequencies of panel modes i and j and of a chosen base or 
reference mode, respectively 


o)-" 



flexibility and stiffness eigenvalues, respectively (eqs. (14) and (15)) 
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ANALYSIS 


Statement of the Problem 

Figxire 1 shows the yawed panel and its coordinate system. A minor difference 
from the analysis in reference 8 is that a right-hand coordinate system is used, so that 
here the z-axis is positive up. Correspondingly, deflections, downwashes, slopes, and 
forces have positive senses upward, but no signs in any of the equations need to be 
chained. 

The usual assumptions are made of linear force -deflection relations, flow over only 
the upper surface, and no aerodynamic effects induced bfetween top and bottom panel sur- 
faces. No aerodynamic perturbation effects on the bottom surface of the panel are 
considered. 

The governing differential equation of motion is that for a uniform isotropic panel 
with no in-plane loading. Two types of structural damping are accounted for: (1) the 
solid-friction mechanical hysteresis type, and (2) the viscous-structural type. The equa- 
tion of equilibrium of forces is given as 

DV^z + iggf) + 0 + Gg II - Ap(x,y, r) = 0 (1) 

where the biharmonic operator is 

s -L ^ + 2 + J:. ^ (2) 

£4 3x4 ^2^2 3x23^2 ^4 

and where D is the panel flexural stiffness, m^ is the panel mass per xmit surface 
area, z is the deflection, r is time, gg£ is a structural damping coefficient of the 
solid-friction (sf) type, Gg represents uniformly distributed panel viscous damping as 
introduced in reference 6, and Ap is the net perturbation pressure (positive with z) 
arising from the deflection and motion of the panel. 

As explained in reference 6 the viscous damping constant Gg is expressed for 
any mode i as 


Gg = gi,v‘"i™A 


(3) 


where ojj is the natural frequency of mode i and g. is a structural damping coef- 
ficient of the viscous (v) type. Furthermore, for Gg constant and independent of fre- 
quency, equation (3) implies that for any two modes i and j 
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( 4 ) 


where subscript B denotes any chosen base or reference mode. 

Direct solution of equation (1) is not feasible for the type of panel deflections and 
aerodynamic forces considered. Therefore, the panel deflection at the flutter condition 
is considered to be sufficiently well approximated by a linear combination of selected 
deflection modes. For simple harmonic motion the time and space variations are sepa- 
rated in the usual way 

z(x,y,T) = z(x,y)ei‘^'T 

where o) is the circular frequency of motion, the shape distribution z(x,y) can be 
complex, and i = is the unit of imaginaries, not to be confused with the modal 

index i. This leads to 

z(x,y,T) Hj(x,y,r) 

j 

Hj(x,y,r) = qj(T)hj(x,y) = qje^^'rhj(x,y) (5) 

where qj is the complex amplitude of the generalized coordinate of motion. When the 
deflection mode shapes hj are known that satisfy both the appropriate boundary condi- 
tions and the differential equation (1), with aerodynamic and damping terms omitted, then 


DV 




= -m. 


tSl. 


dr^ 


^ / , 2v, _io)T 


( 6 ) 


where on the right-hand side co takes on its eigenvalue Wj, the normal-mode frequency 

associated with hj. Thus, the elastic-restoring term Dv4hj is representable in terms 
of ojj. Substitution of equation (5) into equation (1) and use of equations (3), (4), and (6) 
lead to the result 


1 

j 





CO 



- Ap(x,y,T) = 0 


( 7 ) 
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where ggj has taken on its modal values gj g£. The frequency w and the coordinates 
q. are unknown, and the perturbation pressure remains to be specified. 

J 

In keeping with the modal analysis the pressure Ap is also expressed in terms of 
a modal series; that is, 

Ap(x,y,T) “ ^ ^ APj(x,y) = ^ ^Pj(x,y) (8) 

i j 

where each Apj(x,y) is a complex distribution function. 

The Galerkin method is chosen to form the flutter-stability equations. Accordingly, 
the terms of equation (7) are multiplied by a mode shape hj^(x,y) and the resulting prod- 
ucts are integrated over the area of the panel. When the modal index i is made suc- 
cessively 1, 2, 3,. . . for all hj^ employed, the result is a system of equations expressing 
equilibrium of energy for a condition of minimum total energy of the panel. For modes 
that are orthogonal with respect to the mass distribution m^, 

m^h.hj dx dy = 0 (i j) (9) 

S 

each of the set of equilibrium equations is of the form 



From the last term of equation (10) 

£w Apjhj dx dy 3 (11) 

S 

which is the generalized aerodynamic-force (element) term. The procedures for 
obtaining are described in the appendix. 
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Solution of the Flutter Determinant 

For the numerical solutions herein, it has been assumed that the individual modes 
used in the analysis are orthogonal with respect to the mass distribution (no inertial 
coupling between modes), and that there is also no significant stiffness coupling between 
modes. Provision is made for investigating two types of structural damping: (1) solid- 
friction, structural-hysteresis damping characterized by a coefficient gj g^, and 
(2) viscous-structural damping characterized by a coefficient g. . From equations (10) 

J y » 

and (11) the set of equations that express the dynamic equilibrium of motion is given as 



L 




where 


Mj = fw 





dx dy 


(i = 1,2,. . .) (12a) 


(12b) 


and where g^ gj, as introduced in reference 8, is considered as made up of two parts; 
that is. 


^i,sf ®i,sf''‘^sf 


(12c) 


where g^ g^ remains the assigned or measxired value for mode i and gg^ is a modal- 
independent mathematical convenience that is an aid in interpreting eigensolutions. 

Preliminary to a solution of equations (12a), they are each divided by and a 
substitution for Q- is made. The result is 

where 

Qlj = Qij/(pv2jw) 

Two types of solutions of the flutter equations, one adapted to each of the two types of 
structural damping, are described in the following discussion. 


."b 

W B,V 




(i = 1,2,. . .) (12d) 
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The flexibility solution.- A type of solution like that in reference 8, that is suited to 
the presence of solid-friction structural damping only, is termed here the "flexibility 
solution" because the dominant eigenvalue is that of the most flexible mode - the one 
with the lowest of the frequencies cup To obtain this solution based on equations (12) an 
asymptotic expression for 1 + ig^ 


^ + ^i,sf = 1 + %sf + ~ (1 + (l + iggf) 


(13) 


is substituted, g„ is set equal to zero, and equation (12d) is multiplied by 

J3,V 

<1 = 1 . 2 , •■•) ( 11 ») 

where cug is a chosen base or reference frequency, the flexibility eigenvalue is 



(14b) 


1 ^ pf^w ^ . 

M m^£w m^ 


Mq = m^£w 


(14d) 


and k s cu£/V is the reduced frequency. Note that g^ can be different for each 
mode i. 

The stiffness solution .- A type of solution that is suited to analysis of the effects of 
viscous-structural damping, and can simultaneously accoimt for the presence of solid- 
friction damping, is termed here the "stiffness solution" because the dominant eigenvalue 
is that of the stiff est mode - the one with the highest of the frequencies cu^. To obtain 
this solution the modal-independent coefficient gg^ of equation (13) is dropped and equa- 
tion (12d) is divided by cog^. Thus, 



(i = 1,2,. . .) (15a) 
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obtained where the stiffness eigenvalue J2g is 


SI 


s 




(15b) 


As in the flexibility solution gj can be different for each mode i. The eigenvalues 
appear only on the principal diagonal and are readily calculated by an eigenvalue computer 
subroutine. 

Solution techniques .- With aerodynamic forces derived from unsteady potential flow, 
flutter boundaries cannot be calculated directly or explicitly. A search technique must 
be employed, and this is illustrated for a case with the following panel and flow parameters: 
simply supported edges, f/w = A = 10°, two assumed natural beam modes in both 

length and width directions, gf = Oj M = 3.0, and k = 0.25. Equations (14) for the 
flexibility solution and equations (15) for the stiffness solution are solved by using standard 
complex-eigenvalue computer routines, since the eigenvalues 0^ and fig appear only 
on the main diagonal. 

For an initial assumption of 1/ji = 0 the eigenvalues are those for the assumed 
modes in a vacuum. The results are the points at the "foot" of each curve (i.e., at X = 0 
and 1/jn = 0) in all the parts of figure 2(a) from the flexibility solution and of 2(b) from 
the stiffness solution. The labeled values of p,q are the number of half-waves in each 
assumed mode lengthwise (p), and widthwise (q). As 1/ju is monotonically increased 
the traces that are plotted are w* and gg^ ^or gg against X on the left, and stiff- 
ness parameter and gg£ ^or gg y) against 1/jU on the right. The frequency ratio 
coy/ojg and gg^ ^or gg are obtained from the eigenvalues. The dimensionless fre- 
quency parameter co* ^and cOg* and are defined by 


r' 



tjj* 


CO 

f ® 

III 

o -* 


Cl).* 

1 

I 

c^i 

V ^ 




(see appendix A of ref. 8) and oo* is calculated by 


(jO* = 0 )„* — 

B cOg 


(17) 
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The stiffness parameter is 


V ~ 


( 18 ) 


For the stiffness solution, 1/p is obtained from the assumed product (l/p)^a)^o)gj 
and the calculated The dynamic -pressure parameter is calculated according to 

equation (B2) of reference 8 



where j is a function of panel -edge support and 
edges it becomes 


(19) 


£/w, and for simply supported 



( 20 ) 


A sufficient increase of 1/p causes g^j and gg y to pass from negative 
(indicating a stable motion) through zero (neutrally stable, or borderline flutter point) to 
positive values (unstable). The points with ggj = 0 in f^ure 2(a) and with gg y = 0 in 
figure 2(b) are common in all respects, but they have been reached in two different planes. 
Figure 2(a) is in the plane gg y = 0, whereas figure 2(b) is in the plane gi f = 0. (The 
latter could have been in planes equal to an array of constants - possibly a differ- 

ent constant for each mode i - had the analyst so desired.) 

There is a solution difficulty to be overcome at the static aerodynamic limit, k = 0, 
of the imsteady aerodynamic theory. In the equilibrium equations the ratio (l/ju)/k2 
appears. This ratio remains finite when both k and 1/fi approach zero. In order to 
use the same solution framework for calciQating a flutter value for X, an arbitrary non- 
zero value is assigned to k in the equilibrium equations while using the Q.j for k = 0. 

RESULTS AND DISCUSSION 


Introductory Remarks 

The effects of solid-friction and viscous structural damping are illustrated for the 
same panel characteristics studied above in the section ’’Solution Techniques,” namely, 
simply supported edges, i/w = \/o73, A = 10^, and two assumed modes in both length and 
width directions. Effects of solid friction alone are obtained from the flexibility solution. 
Effects of viscous friction and of combinations of viscous and solid-friction damping are 
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from the stiffness solution. Comparisons are shown for k = 0 representing the static- 
aerodynamic approximation, for k = 0.01 representing the quasi-static aerodynamic 
approximation, and for k = 0.25 representing the use of the complete unsteady aero- 
dynamics. The majority of the studies are done with M = 3.0, as typical of the Mach 
number range where coupled-mode flutter is encountered, but one study is reported for 
M = 1.3 where "single-degree-of -freedom'* flutter can be encoimtered. 

In the "Solution Techniques" section the search procedure for finding flutter boimd- 
ary points was described with plots of OJ* and g^j ^or gg against X and of stiff- 
ness parameter and gg£ ^or gg against l/p. The same types of plots are used in 
presenting the results that follow. 

Comparison of Damping Effects Using the Static Aerodynamic Limit (k = 0) 

It was pointed out at the end of the "Solution Techniques" section that, to obtain 
restilts for k = 0, the are calculated for k = 0, but in the equilibrium equations k 

is assigned a nonzero value. Although any value can be assigned, k = 0.01 was assigned 
in the present work. This choice necessarily affects the range of 1/p that is selected 
during the search. 

Solid-friction effects from a flexibility solution .- The parameter traces are shown 
in figure 3(a). An initial choice of 1/p = 0 gives the natural-mode -in-a-vacuum points 
at the foot of each curve where the p,q values are labeled for reference. As 1/p is 
gradually increased the parameter traces move as indicated by the arrows. A small 
value of 1/p results in the four points indicated by the circles. The four circles fall 
at a constant 1/p on the plot of stiffness parameter against 1/p; but on the plot of oj* 

n 

against X the four circles fall on a parabola, and X is a constant multiplied by {(*}*). 
For a higher chosen 1/p the four points indicated by the squares resvilt and they fall 
on a lower parabola on the plot of w* against X. This 1/p was chosen because two 
frequencies have coalesced at an instability point with X = 331. (Two of the squares ' 
coincide at the frequency coalescence of the p,q =2,1 and 2,2 modes.) A still higher 
value of 1/p gives the diamond points. The triangle points result from a yet higher 
value of 1/p that was selected for its closeness to the frequency coalescence of the 
p,q = 1,1 and 1,2 modes with X = 168. Note that the parabola through that coalescenCe- 
point pair is tangent to the curve of against X. 

With static aerodynamics the only way a flutter instability can be reached is by a 
frequency coalescence, and this occurs for the present case with X = 168, which is the 
flutter value. 

Viscous-friction result from a stiffness solution.- The parameter traces are shown 
in figure 3(b) and the natural-mode values of p,q are labeled at the foot of each curve. 
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A small chosen value of (l//i)^aycOg^ gives the parameter values Indicated by the 
circles, and certain successively higher values give the squares, the diamonds, and the 
trinities. Note that each set of symbols falls on a vertical line on o)* plotted against X, 
but the syrhbols fall on a parabola on stiffness parameter plotted £^ainst 1/p. Compar- 
ison with the flexibility- solution resvilts in figure 3(a) reveals the contrastii^ behavior of 
any set of points (such as the circles). 

On the plot of w* against X the lower frequency pair of squares falls just barely 
above the flutter associated with frequency coalescence at X = 270 of the modes that 
originate at p,q = 1,1 and 1,2, and this frequency coalescence occurs at a point of tan- 
gency with the vertical line X = 270. The diamond points fall at an intermediate value 
of X. At a higher X a pair of the triangle points falls very nearly at the frequency 
coalescence of the modes originating at p,q = 2,2 and 2,1. 

It is significant to observe the two different results for k = 0. From the flexibility 
solution with gg y = 0 with ^gj + g^j^ variable including — 0, flutter occurs at a 
frequency coalescence at which X = 168. (See fig. 3(a).) From the stiffness solution 
with g. gj = ggj = 0 and with gg ^ variable including gg y frequency-coalescence 
flutter occurs at a frequency coalescence with X = 270. This finding is discussed further 
subsequently. 

Viscous-plus-solid-friction result from stiffness solution .- In order to study com- 
bined effects of viscous and solid friction, especially for small amouts of solid-friction 
damping, gj^ was set equal to 0.01 (all modes). The results from a stiffness solution 
are shown in figure 3(c). At the crossing point g„ „ = 0, the same result is obtained 
(X = 168) as from the flexibility solution with v ^ ® ^sf ~ shown in fig- 

ure 3(a). Small increments of g^ ^ above 0 are strongly stabilizing. For gg ^ = 0.001, 
X = 176; and for „ = 0.01, X = ^29. ^ 

r>,V 

A separate solution for g^^ = 0.001 was made in order to approach more closely 

the limit gj gf 0- The flutter X was again 168, and increase of gg ^ was even more 
strongly stabilizing than in the presence of gj gf = 0.01. 

Although such solutions that are very close to gj gf = 0 do not constitute a formal 
mathematical proof of the limiting value of X, it is strongly believed that X = 168 is the 
limit value as g^ gf 0 within the stiffness solution because of the smooth and regular 
behavior of all the forces acting. 

Summary plot of effects of solid -friction and viscous damping.- Figure 4 represents 
a three-dimensional plot of X against both gg ^ and ^g^^ The curve seg- 

ments in the three planes gg y = 0, gj gf = 0, and g^ = 0.01 are from figures 3(a), 
3(b), and 3(c), respectively. The discontinuous result for X for g^^ — 0 compared 
with gg£ s 0 is shown. 
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Comparison of Damping Effects Using Quasi-Static 
Aerodynamics (k = 0.01) 

A very low value of k (0.01) is used in order to study the damping effects in the 
vicinity of the limit k — 0, where the quasi-static aerodynamic approximation is most 
valid. 

Solid-friction e ffects fr om a flexibility solution. - The parameter traces are shown 
in figure 5(a) and, being from the flexibility solution, are for the condition (and in the 
plane) gg ^ = 0. In order to display better the pertinent trace behavior, the range of the 
ggf scale is between ±0.02, which is much less than in figxares 2 and 3. The points on 
the traces at certain successively higher 1/fj, values indicated by the circles, squares, 
diamonds, and triangles are at the same 1/jix values as in figm-e 3(a). The trace or 
mode that reaches the flutter point is the one originating at p,q = 1,2, and for g^^ > 0 
the curves on all figure parts are a series of dots for the purpose of identification. At 
the crossing point ggj = 0, \ = 270. For g^^ = 0.01, X drops sharply to 215 indi- 
cating the strong destabilizing effect of solid friction in this parameter range. (The 
other trace that reaches an instability at a higher X originates at p,q = 2,2, and for 
ggf >0 it is identified by a row of x's on each figure part.) 

Viscous -friction resul t from stiffness solution .- The parameter traces for gj gf = 0 

are shown in figxme 5(b). The rai^e of gg ^ is restored to ±0.4 in order to display 

better the trace behavior. At the crossing point gg y = 0, X = 270. For a small or even 

substantial increase in g„ „ the increase in X is insignificantly small. 

is,v 

Viscous-plus-solid-friction re sult from a stiffness solution .- Values of gj^ for 
all indexes i are set equal to 0.01, the same as for k = 0 described previously, and 
the results are shown in figure 5(c). At the crossing point gg ^ = 0, X = 215 which is 
the same as for g^^ = 0.01 in figure 5(a), as it should be, since they are a common point 
reached in two different intersecting solution surfaces. An increase of gg ^ above 0 
is stabilizing: for gg y = 0.01, X = 236. 

Summary plot of effects of so lid-fric tion and viscous dampi ng.- Figure 6 represents 
a three-dimensional plot of X against both gg ^ and ^g^ + g^^^. In comparison to 

figure 4 with k = 0, here the discontinuous behavior has disappeared and there is a steep 
negative slope with varying ^g- The curves in the three planes gg ^ = 0, 

gj = 0, and gj^ gj = 0.01 are from figures 5(a), 5(b), and 5(c), respectively. 

Although a formal mathematical proof is lacking, it is concluded that there is 
a discontinuous result for x from the flexibility solution for k — 0; namely, for 
gi g£ + gg£ = 0, X = 168 for k = 0, but X = 270 for k — 0. This conclusion is reached 
because of the smooth and regtilar behavior of the aerodynamic forces in the vicinity of 



k — 0 and from observation of how an Instability (flutter) point is reached in the solution 
processes. 


Comparison of Damping Effects Using 
Complete Unsteady Aerodynamics 

A result of using the complete unsteady aerodynamic theory is shown in figures 2(a) 
and 2(b) for k = 0.25. This figure was used earlier to illustrate the technique of solu- 
tion. The flutter point with \ = 293 is at the crossing point g^j = 0 in f^ure 2(a) and 
at the crossing point gg y = 0 in figure 2(b). An increase of g^£ to 0.02 predicts a 
small decrease of X to ^286. An increase of gg ^ to 0.02 predicts a small increase 
of X to 302. 

At the common point, gg£ = gg y = 0, X = 293, l/p = 0.0721, = 0.1196, 

and- the product (l/p,)(o)j£^V) = 0.00’862. Figure 7 shows a flutter boundary in the 
plane 1/jU plotted against stiffness parameter as used in reference 3 and in 

subsequent publications. The solution points for k = 0.25 and 0.01 discussed previously 
are indicated. A variable thickness of panel with a fixed panel material, altitude, and 
Mach number determines a hyperbola (l/jLt)^a>j^£^V) = Constant. The hyperbola of short 
dashes is for aluminum at standard sea level and with M = 3, for which (l//i)(o)j^£^vj 
= 0.00851. Increasing thickness is to the right toward the stable region. The intersec- 
tion of the hyperbola with the flutter boimdary determines the thickness ratio required to 
prevent flutter. 

For panel and flow parameters that result in a flutter boundary characterized by a 
slow variation of X and of with varying k and altitude, extrapolations and 

interpolations to obtain additional points on a flutter boundary can be obtained on the basis 
of assumed constant X and w/wg. The associated proportionalities are 


1/M _ (Vm)(c0i£/V) ^ (1/m)KVV) 


( 21 ) 


from which 


k oc 




( 22 ) 


When an analyst has tried a value of k and obtained a point on a flutter boundary that 
does not fall on a desired hyperbola of panel material and altitude, he can estimate his 
next value of k from equation (22). 
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Comparison of Damping Effects at M = 1.3 

All of the results discussed previously are for M = 3.0, as representative of a 
Mach nimiber well above any transonic or low supersonic effects. In this higher range 
of M, panel flutter is characterized by strong coupling of the assximed deflection modes 
used in a modal analysis. In the low supersonic range, for M less than about \J2, and 
for panels with £/w less than about 2 to 4, a different type of flutter is predicted. The 
flutter motion is dominated by a single mode and is often termed "single -degree -of - 
freedom flutter," although the other modes are present in small proportions. 

Results calculated for the same panel discussed previously, but with M = 1.3, are 
presented in figure 8. As indicated by the key and the labels, the solid boundary on the 
right is the flutter boundary for gg y = ggf = 0- For all the lower part of that boundary 
the flutter motion is predominantly the p,q = 2,1 mode. An increase of viscous damping 
alone to g^j „ = 0.01 has a significant stabilizing effect by shrinking the flutter region to 
within the bovindary of short dashes. An increase of solid-friction damping alone has a 
comparatively greater stabilizing effect by shrinking the flutter region to within the 
boundary of long and short dashes. The solid, nearly parabolic -shaped boundary to the 
left is an instability boundary characterized by strong modal coupling, and it is little 
affected by structural damping except close to the origin where k is very small. A 
sufficient increase of gg ^ and/or g^^ shrinks the flutter region that is dominated by 
the single -degree -of -freedom (p,q = 2,1) mode, imtil the coupled-mode instability bovind- 
ary becomes the flutter boundary. 

Comments on Convergence With Respect to the 
Number of Modes in the Modal Analysis 

All of the results of the present report were determined by using four modes in the 
Galerkin modal analysis - two modes lengthwise and two widthwise even though these 
results are not particularly well converged. The use of four modes in both directions 
gives a X roughly 30 percent higher, and it was established that this four-by -four-mode 
result is well converged, at least for A = 0° and 90°. The effects of solid-friction and 
viscous structural damping are the primary subject of this report, and these effects are 
the same from the two-by-two-mode results and from the four -by -four-mode results. 

The 16 -mode results are so much more complicated to present and describe that the 
decision was made to present the 4-mode results. 

Review of Results 

The two types of structural damping, solid friction and viscous, are shown to have 
mostly contrasting effects on panel -flutter boundaries for the panel and flow parameters 
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investigated. Solid-friction damping is strongly to moderately destabilizing, whereas vis^ 
cous damping is stabilizing or ineffective. This stabilizing effect of viscous damping on 
panels contrasts with the findings of Ellen in reference 7 for two-dimensional membranes 
based on quasi-static aerodsmamics. Figure 1(b) of reference 7 shows a marked destabi- 
liy. ing effect of both viscous and solid-friction (i.e., hysteretic) dampir^. Since the effects 
of the two types of structural damping are foimd to be so different herein, there cannot be 
an ’’equivalent viscous damping” for solid -friction structural damping in panel flutter. 

Flutter solutions were obtained for a particular yawed panel for M = 3.0. With use 
of the static -aerodynamic approximation (k = 0), the stiffness solution gives a discontin- 
uous result for \ in the vicinity of Sj gf “ 0 the plane gg y = 0- For g. = 0, 

X = 270; but for g^ gf ^ = 168. Another discontinuous result is obtained from the 

flexibility solution in the vicinity k » 0 for g^^ = gg v ^ k s 0 (i.e., static 

aerodynamics), X = 168; but for k — 0 (i.e., quasi-static aerodynamics), X = 270. 

The complete unsteady aerodynamic forces were used for k = 0.25 in order to 
calculate a point on the flutter boundary that corresponds to aluminum panels at sea level. 
All of the flutter boundary results for M = 3 are characterized by a flutter motion of 
strongly coupled assumed deflection modes. 

The low Mach number region was studied for the same panel with M = 1.3. Here, 
the flutter boundary is characterized for low structural damping by a flutter motion that 
is dominated by just one of the assumed modes (”single-degree-of -freedom flutter” 
herein termed ”DOF”). This single DOF boundary is not at all characterized by a nearly 
constant X. Addition of both solid-friction and viscous damping is strongly stabilizing, 
and a sufficient increase of either type of damping rapidly shrinks the single DOF flutter 
region until a coupled-mode flutter boundary is encountered. This latter boundary is 
like the coupled-mode bovindary for M = 3 in that X varies slowly with 1/p, and the 
boundary is little affected by further increases in damping, except at very low l/p. 

CONCLUSIONS 

An equilibrium equation has been given that accounts for both solid-friction and 
• viscous structural damping in panel -flutter analysis. A solution technique has been 

developed for routinely calculating the effects of both types separately and in combination. 
^ The solution technique is usable with aerodynamic forces from general unsteady potential 

flow. A Galerkin modal analysis of a rectangular yawed panel at two Mach numbers M 
of 3.0 and 1.3 leads to the following conclusions; 

1. The additions of solid -friction and viscous structural damping to the panel have 
opposing effects on the coupled-mode flutter boundary: solid-friction damping is mostly 
destabilizing, but viscous damping is stabilizing or ineffective in every case investigated. 


17 



I 


II I II n I II 1 1 


1 


The relative effects of structural damping are the greatest at the lowest mass ratios 
(air mass to panel mass) and become smaller as mass ratio increases. 

2. Since the two types of dampii^ result in opposing effects at M = 3.0, there can- 
not be, in general, an "equivalent viscous" damping to solid-friction damping. 

3. Two discontinuous results for flutter dynamic -pressure parameter are obtained 
at M = 3.0 in the vicinity of reduced frequency k « 0 and in the vicinity of solid- 
friction damping g^j » 0. 

4. At a low supersonic Mach number, M = 1.3, the flutter boundary is a "single- 
degree -of -freedom" flutter boundary for the important rai^e of mass ratios and for low 
structural damping. Both types of damping are strongly stabilizing for the single -degree - 
of -freedom flutter. A sufficient increase of structural damping shrinks the single- 
degree-of -freedom flutter region until a coupled-mode flutter boundary is encovintered, 
for which the structural damping effects are similar to those for M = 3.0. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
August 2, 1976 
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APPENDIX 


GENERALIZED AERODYNAMIC-FORCE ELEMENTS Q^j 
The generalized aerodynamic -force element is, as given by equation (11), 

Qjj = £w jy Apjhj dx dy (Al) 

S 

The lifting-pressure difference Apj is that from the modal series of equation ( 8 ) 

Ap(x,y,r) » ^ ^ Apj(x,y) (A2) 

j 

The lifting pressure is obtained from the perturbation velocity potential <p(x,y,T) by 


Ap(x,y,r) = -p (A3) 

By using a modal series for <p that is consistent with the series for Ap and for 
deflection z, 


q * 


and 


Ap 


{x,y,t) = -pVe^^^'’’ Y -/(cos A — 7 ^ + sin A i 

a(£x) a(wy) 




which gives for the individual Apj 


a(fx) 


d(p. 


Ap. = -pV (cos A — 7 - 4 t + sin A . J . + <Pi 1 

a(wy) V ij 


(A4) 


(A5) 


(A 6 ) 
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where 


?>j(x(x,y).y(x,y))=j^j3’ 


Wi(«,,)e-‘“‘=‘-»cos(|R) 

V r" 


dri 


(A7) 


wherein; 

S* region of integration, the panel area within the forward-facing Mach charac- 

teristic cone with apex at x,y 


x,y streamwise and cross-stream dimensional coordinates, respectively, with 

same origin as x,y: 

X = fx cos A + wy sin A 

y = wy cos A - £x sin A 


^, 7 ] dummy variables for x,y, respectively 


R 

(5 


= \/(x - - /s2(y _ rjf 

= 


and the downwash amplitude is 


8h, 


9h, 


ct)£ 


w. = Vicos A — i + — sin A — ^ + i 
J \ 8 x w ay V 

Substitution of equation (A6) into equation (Al) gives 



(A8) 
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Integration by parts of the terms that involve sin A and cos A yields 



(A9) 


Where panel edges can displace vertically, the terms hj<^j evaluated at those edges can 
be significant; but for the usual situation of hj^ = 0 at the edges, Q.j is determined by 
the remaining terms 


Qjj = pVw 


= pVw 



8h, 


9h, 


cos A + 4 sin A — ^ - i ^ hs ] <p. dx dy 

8x w - ay V V J 



(p^ dx dy 


(AlO) 


where the velocity potential (p^ is that of equation (A7) and w-’yv is the complex con- 
jugate of the downwash ratio (Compare with Wj of eq. (A7).) The elements Qy 

as given in equation (AlO) are used in the present report. 
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Dynamic-pressure parameter, X 


Mass ratio, 1/m 




(a) Flexibility solution. 

Figure 2.- Parameter traces from flutter solution for panel with £/w = \/0.3, 

A = 10°, M = 3.0, and k = 0.25. 
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Figure 3. 
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(a) Flexibility solution. gf = 0- 

Parameter traces from flutter solution for panel with £/w = /o.3, 
A = 10°, M = 3.0, and k = 0. 
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(b) Stiffness solution, g. = 0. 
Figure 3.- Continued. 
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(c) Stiffness solution, g. = 0.01. 
Figure 3.- Concluded. 



Figure 4.- Dynamic -pressure parameter X against g^^ with gg y = 0, 
and against gg ^ with g^j = 0 for panel with £/w = \/0.3, A = 10°, 
M = 3.0, and k = 0. 
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Dynamic-pressure parameter, X Mass ratio, l/p 


(a) Flexibility solution, gj gf = 0- 

Figure 5.- Parameter traces from flutter solution for panel with f/w = \Zo.3, 

A =10®, M= 3.0, and k = 0.01. 




Dynamic-pressure parameter, X 


Mass ratio, 1/fj. 



Dynamic-pressure parameter, X Mass ratio, 1//li 


(c) Stiffness solution, = 0.01. 

Figure 5.- Concluded. 
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Figure 7.- Flutter boundary. Mass ratio 1/ji against stiffness 
parameter for panel with £/w = \A).3, A = 10 , 

g . = 0, and M = 3. 
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